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Abstract 

The many-particle spectrum of an isotropic frequency gap medium doped 
with impurity resonance àtoms is studied using the Bethe ansatz technique. 
The spectrum is shown to contain pairs of quantum correlated "gap excita- 
tions" and their heavy bound complexes ( "gap solitons" ) , enabling the prop- 
agation of quantum information within the classically forbidden gap. In addi- 
tion, multi-particle localization of the radiation and the medium polarization 
occurs when such a gap soliton is pinned to the impurity atom. 
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Light localization is a classical effect predicted Q to occur in strongly scattering dielectric 
microstructures. In the context of photonic band gap (PBG) materials nonclassical 
forms of localization such as photon-atom bound states have been predicted [ffl when the 
resonant transition frequency of an impurity atom lies within a gap. This bound state is an 
eigenstate of the quantum electrodynamic Hamiltonian for a realistic PBG crystal exhibiting 
a general anisotropic photon dispersion relation. In this state, a virtually emitted photon 
may tunnel many wavelengths away from the atom before being reabsorbed, leading to non- 
Markov memory effects || in collective light emission from many àtoms. It was recently 
shown that an effective model || describing both isotropic PBG systems and frequency 
dispersive media (DM) |7|] doped with resonance àtoms exhibits hidden integrability and 
is diagonalized exactly f||| by means of the Bethe ansatz technique [f|. This suggests the 
possibility of a rich multi-particle spectrum in real physical systems exhibiting a frequency 
gap, when such systems are doped with impurity àtoms. 

In this paper, we demonstrate the existence of nonclassical states of light which may be 
generated, for instance, through the interaction of an external làser field with an impurity 
atom placed within a polariton gap [|K]] of a DM. In addition to ordinary polaritons and their 
bound complexes ("ordinary solitons") occurring outside of the gap, the sub-gap spectrum of 
the system is shown to contain propagating pairs of correlated "gap excitations" and their 
heavy bound complexes ("gap solitons"). The individual gap excitations comprising the 
pair are correlated such that the probability amplitude of finding them far apart decreases 
exponentially with the ratio of their separation distance to the classical penetration length 
of the radiation into the medium. In addition to heavy gap solitons propagating within 
the gap, the spectrum contains multi-photon localized states pinned to the atom. Under 
external perturbations a pinned gap soliton may dissociate into propagating gap excitations. 
We evaluate the dispersion relations, the effective masses, and the dissociation energies of 
quantum gap solitons and we show that they are stable with respect to weak perturbations. 
Our results demonstrate a clear distinction, at the quantum level, between fermionic gap 
systems (such as electrònic semiconductors) and bosonic gap systems. In a semiconductor, 
propagation within the energy gap is strictly forbidden. In the bosonic gap, however, certain 
nonclassical many-body gap states are allowed to propagate. This may have important 
consequences for the transmission of quantum information within a bosonic gap medium. 
Although the Bethe ansatz method made use of the isotropic one-particle dispersion relation, 
qualitative similar results may hold in a realistic, anisotropic PBG material which also 
exhibits a small nonresonant Kerr nonlinearity [JTT . 



Consider the model Hamiltonian H = H + V, where 

H = cj 12 (a z + 1/2) + / ^up*(u)p(u) (la) 
y = -y/ï í %{u)\p{u)o-++p\u)o-]. (lb) 

JC Z7T 

Here, uj\2 is the transition frequency of the two-level impurity atom placed within a bosonic 
frequency gap medium. The spin operators a = (a x ,a y ,a z ), (j ± = a x ± io v satisfy the 
Standard commutation àlgebra [er\ er J ] = €ijkO~ k and act on the atòmic variables of the 
system. The operators p\uj) (p(u>)) create (annihilate) bosons of frequency ui in a specific 
(elèctric dipole) spherical harmònic state and satisfy the àlgebra [p(o>),pt(ü/)] = 2·n5{uj — ui'). 
The integration contour consists of two parts, C = C_ © C + , where C_ = (0, Q±) and 
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C + = (0||,oo) correspond to the lower and upper branches of the medium excitations, 
respectively. The interaction term (lb) describes emission and absorption of bosons by the 
atom. Here 7 = 4ujf 2 d 2 /3 is the inverse lifetime of the excited atom in free space with 
the dipole transition moment d, while the atòmic form factor z(u) contains the information 
about the polariton spectrum. Since the polariton gap persists for small wave vectors, k — > 0, 
we can neglect anisotropies of the underlying iònic crystal and take the dispersion relation 
to be isotropic. Also, we choose units in which h — c—1. 

In the one-particle sector of the full Hilbert space, it is straightforward to verify p|,|6|,|T2 



that when a> 12 lies inside the frequency gap, one of the eigenstates of the Hamiltonian H 
describes a polariton-atom bound state. The multi-particle sector of the Hilbert space can 
be studied by the Bethe ansatz technique. The Schròdinger equation (H — E)\^n) = is 
solved exactly ||[8] due to the two-polariton factorization of the multi-polariton scattering. 
Imposing the periòdic boundary conditions (PBC) on the iV-polariton wave function leads 
to the following set of Bethe ansatz equations (BAE): 



e 



ik , hj - 0/2 _ £L hj -h-W 



hj + if3/2 M hj — hi + if3 



n J : p) 



which completely determine the iV-particle spectrum of the model (1). Here, E = YljLi is 
the eigenenergy and Uj are polariton freqüències which solve Eqs. (2). Also L is the radius 
of a sphere centered at the atom on which we apply PBC and then take the limit L — > 00. 
In the case of a DM, the polariton momenta, kj = k(ujj), and "rapidities" , hj = h(uüj), are 
expressed as, 

k(u) = um(u), h{u) = (—) - ^ 12 , (3) 



with the refractive index n(ui) = \Js{u) and the dielectric permeability of the medium 
e(uú) = (u) 2 — íí|)/ {uj 2 — ííj_). Here Ü12 — UJ12 is the Lamb shifted atòmic transition frequency. 
The parameter j3 = 7/^12 appears in both the polariton-atom scattering (l.h.s. of Eqs. (2)) 
and in the effective polariton-polariton coupling (r.h.s. of Eqs. (2)) caused by the polariton- 
atom scattering. 

As L —>■ 00, apart from real solutions, Eqs. (2) admit complex ones, in which the 
rapidities hj are grouped into the Bethe "strings" . In this paper, we confine ourselves to the 
case when all iV rapidities are grouped into a single string hj = H + i%(N + 1 — 2j), j — 
1, . . . , N with a common real part ("carrying" rapidity) H. This is a solution of BAE if and 
only if the imaginary parts of rapidities hj and corresponding momenta kj have the same 
sign: 

sgn(Im hj) = sgn(Im kj), j — 1, . . . , N . (4) 

This restricts possible magnitudes of polariton freqüències ujj corresponding to the string 
rapidities. It is easy to understand that the necessary condition (NC) (4) determines the 
frequency intervals, in which the effective polariton-polariton coupling is attractive leading 
to bound many-particle complexes (quantum solitons). 

In empty space, an effective photon-photon coupling is attractive for all freqüències of 
physical interest, and a Bethe string in the space of rapidities is mapped to a quantum 
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soliton in the space of freqüències |Ï3| , ujj = Q + i^(N + 1 — 2j) with E = QN. This also has 



a string structure. To avoid a possible confusion in what follows we use the term "string" 
for solutions of BAE in the /i-space and the term "soliton" to refer to string's images in the 
lo- and fc-spaces. 

In the medium, analytical continuations of the functions k{uS) and h(u>) in the complex 
u>-plane depend essentially on the position of the real part of the frequency with respect to 
the medium gap. We start with the case when the real part of u lies outside the gap. Let 
uj = X + ir] and the real part X E C. For 77 A, the functions k(u) and h(u) are then 
represented as k(u) = k(X) + irjk'(X) and h{u) = h(X) + ir)h'(X). Since k'(X) = dk(X)/dX is 
positive for all À, NC leads now to the condition h'(X) = dh(X)/dX > 0. This is met only if 
À E CL. Therefore, the effective coupling is attractive only between polaritons of the lower 
branch, and the soliton image of the Bethe string is given by ujj = Q + ~r(íí)(iV + 1 — 2j) 
and kj = K(Q) + |Q(íí)(JV + 1 - 2j). Here K(Q) = Qn(Q), T(Q) = P/h'(Q), Q(Q) = 
/3k'(Q)/h'(Q), and E = QN is the soliton eigenenergy. The common real part of the polariton 
freqüències is found from the equation h(Q) = H, which has a root lying in CL only if 
H < 0. The soliton obtained is quite similar to a vacuum soliton, and we will use the 
phrase "ordinary soliton" to refer to this solution, despite its inordinate behavior on different 
polariton branches. Polaritons of the upper branch are described by one-particle Bethe 
strings with real positive rapidities and do not form any bound complexes. The results 
obtained are clearly vàlid if uu E C- © G, where G = (Q±,Q\\). If üüu lies above the 
gap, E C + , the effective coupling, including interbranch one, becomes attractive and 
admits both ordinary solitons in each branch and unusual "composite solitons" containing 
polaritons of different branches. 

Now let us look for an image of a Bethe string, provided the real parts of all the freqüèn- 
cies ujj lie inside the gap. Let u = £ + irj and £ E G. To find the analytical continuations 
of the functions k(u) and h(u) to the complex u space, we need first to fix an appropriate 



branch of the function n(u>). Let n(£ ± iO) = ±w/(£), where u(Ç) = y\s(Ç)\. In this case 
k(u) = sgn(77)[-77/í'(£) + ïk(£)], where k(£) = £//(£)• Also h(u) = sgn(r}) [7}<//(g) - i<j)(£)] 
where </>(£) = (£ — ^12) /(0 5 an d /(0 = ^12 [^ 3l/5 (0] -1 - Since the function is positive, 
NC leads to the condition </>(£) < 0. It means that allowed gap excitations exist only for 
£ G CU12). Because of a strong nonradiative relaxation in the medium in the vicinity of 
the frequency fl±, we focus our studies on gap states of physical interest lying in the vicinity 
of the atòmic frequency uu- The remaining analysis is simplified by linearizing the function 
</>(£) at the point £ = 0^12, 0(0 — a(£ — uiu), where a = /(a;^). 

Now we are able to map a Bethe string to corresponding gap excitations. We start 
with the simplest case of a two-particle string, N = 2. Its complex conjugated rapidities, 
h = H + i/3/2 and h* = H — i/3/2, are mapped to the corresponding pairs of the complex 
conjugated freqüències, uj = £ + irj and oj* = £ — irj, where the imaginary part is assumed 
to be positive, i] > 0, and momenta, k = q + i/í(£) and k* = q — i/í(£), where q = —t]k'(^). 
The real and imaginary parts of freqüències are expressed in terms of the string parameters, 
£ = üüi2 — (3 /2a and rj = H/a. In the spherical harmònic formalism introduced previously 
||, the real part of the particle momenta q must be positive. Since is negative, it 

follows that q = T]\k'\. Consequently only a string with a positive carrying rapidity, H > 0, 
is mapped to gap states of physical interest. 

The expressions obtained describe a novel, quantum correlated state of two gap excita- 
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tions. Two gap partides comprising the pair are "confined" to travel together. They do 
not exist separately from each other, unlike polaritons of ordinary solitons. Moreover, the 
confined state cannot be treated as a bound state of two polaritons from different branches 
(like a Wannier-Mott exciton in semiconductors) , because, under the condition ui 2 G G, the 
interbranch polariton-polariton coupling is repulsive. The spatial size of a pair, 5 ~ 
is nothing but the penetration length of the classical radiation field with the frequency 
íü = Ç G G into the medium O. Since the wave function of a single gap particle is unnor- 
malizable, free one-particle gap states in the bosonic gap are forbidden in exactly the same 
way that electrònic propagation is forbidden in a conventional semiconductor gap. However, 
in the case of bosons, the effective particle-particle coupling allows one to construct the 
normalizable wave function of a pair from unnormalizable wave functions of each particle. 
At large interparticle separations, the wave function of a pair has the form 



where the real and imaginary parts of momenta describe, respectively, the motion of the 
center of gravity and the spatial size of a pair. The auxiliary coordinate variables Xj are 
analogous to spatial coordinates along an arbitrary axis passing through the atom. The 
vicinities x < and x > correspond, respectively, to ingoing and outgoing spherical 
harmonies of the polariton field. An angular distribution of the field is determined by the 
specific spherical harmònic polariton state. 

Let us consider now the mapping of a string containing an even number of partides 
("even string"), N = 21, to gap excitations. A pair of complex conjugated rapidities, hj 
and hj, j = 1, . . . , l, is mapped to a pair of freqüències Uj = ^ + irj and u* = ^ — irj 
and corresponding momenta kj = qj + ík(^) and k* = qj — where the real parts 

of freqüències are given by ^ = ujm — (j3/a)(l + 1/2 — j). The upper index, (0), indicates 
that these expressions are derived within the linear approximation for </>(£) in the vicinity 
of £ = Ui2- The expressions obtained describe a bound complex of l pairs of confined gap 
partides with the eigenenergy 



We will use the phrase "gap soliton" to refer to this state of the system. Unlike an ordinary 
soliton, a gap soliton is stable with respect to quite weak perturbations of the system, 
because the dissociation energy, Ud, of a soliton with l pairs into two solitons with l\ and 
h — l — h pairs is positive: Ud = Ej^ + — E^ = 2(f3/a)l\l2- The radial thickness of the 
spherical harmònic soliton pulse is determined by the imaginary part of the momentum k\ 
corresponding to the rapidity h\. Since is monotonically decreasing function, the gap 
soliton size, òf^ ~ ^ _1 (£Í°^), falls with the growth of the number of pairs l. 

Since the effective coupling constant is very small, /3 <C 1, the linear approximation works 
well even for large solitons containing many pairs. But, in this approximation, the soliton 
energy is independent of its momentum. For what follows it is convenient to introduce 
the soliton energy per particle ej^ = Ej ' /2l = co> 12 — (f3/2a)l, and the soliton momentum 
per particle, q = (2/) _1 J2j Qj ~ (v/l)J2j )l- To estimate the first corrections to Eq. 




(5) 




2J2Q° =2cü 12 l-(f3/a)l 2 . 



(6) 
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(6), we have to keep the next term, ï(?7 2 /2)0"(£), of the Taylor series for the function h(uj) 
and the next term in the expansion of the function 0(£) at the point £ = cl> 12 : — 
a(£ — U12) + ò(£ — u; 12 ) 2 , where ò = f'{üJvz) > 0. The freqüències £j are then given by 
£j — £j°' ) = — (b/a)(C,^ — CJ12) 2 + (b/a)rj 2 , while the soliton momentum is still given by g. The 

term, — (ò/a)(£j ' 1 — co> 12 ) 2 , leads to the first order correction to the energy of a motionless 
soliton and determines the width of the soliton band, while (b/a)rj 2 leads to the kinetic 
energy contribution to the total soliton energy in the effective mass approximation: 

e, = e! 0) - A t + q 2 /2m x . (7) 

Here A; = ^(4/ 2 - 1) and m, = ^2 fej Wi^lY are the band half-width and the 
effective mass, which increase with l. At small l, the propagating gap soliton bands are 
very narrow and solitons are very heavy and even motionless at q — 0. But the bandwidth 
increases as l 2 , so that large l solitons are quite mobile when the momentum q becomes 
larger than the range of validity of the effective mass approximation. At arbitrary q, the 
exact equations for the soliton parameters and rjj are given by Reh(£,j,r]j) = H and 
lmh(£j, rjj) = (3(1 — j + 1/2). The solution of these equations requires simple numerical 
calculations. 

Finally, we evaluate the pinning energy of a gap soliton to the atom. In /i-space, pinned 
solitons are described by odd strings with H — > + . The one-particle string, l = 0, with 
H — ► + is clearly mapped to the gap state with u = Ujyi ~ ujyi. This state is nothing but the 
polariton-atom bound state [|4||6l,|ï2"| in the one-particle spectrum of the system. Therefore 
the extra real rapidity of an odd string can be mapped to the gap state with ^ = uji 2 , while 
the remaining complex conjugated pairs of rapidities are mapped to a deformed motionless 
gap soliton with = uJi 2 — ([3/a)(l—j + l). We used here the term "deformed" to emphasize 
that the soliton freqüències now contain the extra term — (3 /2a due to the polariton-atom 
bound state, which deforms the soliton, in contrast to the mobile, even gap soliton whose 
freqüències are given by Since one of the partides of an odd motionless gap soliton is 
bound to the atom, a soliton as a whole is also "pinned" to the atom. Moreover, as H — > + , 
the imaginary parts of freqüències and the real parts of particle momenta vanish. A pinned 
soliton describes a many-particle state of the system, in which the radiation and medium 
polarization are localized in the vicinity of the atom. To evaluate the energy of pinning, 
we need only to compare the energy of a soliton with 21 + 1 partides pinned to the atom, 
= 0012(21 + 1) — (f3/a)l(l + 1), with the sum of the energies of the one-particle bound 
state, Un-, and a motionless gap soliton with l pairs, Ef\ We find that the binding energy 
Ui = —((3/a)l is proportional to the number of pairs. Moreover, the energy required to pull 
a single pair of partides out of a pinned soliton, U\ = (E^/ x + Ef^) — Ei = ((3/a)(2l — 1), 
is even greater than the energy required to pull off all / pairs. Therefore, the state of a 
pinned soliton is stable with respect to quite weak perturbations of the system and its 
stability increases with l. We mention finally that, for H > 0, it is also possible to construct 
composite solitons, consisting of an odd number of bosons, which correspond, physically, to 
a bound state of a gap soliton and a polariton in the upper branch. 

In summary, we have shown that the isotropic dispersive medium doped with an impurity 
atom exhibits a rich many-particle spectrum containing heavy, mobile, gap solitons as well 
as pinned solitons. Mobile gap solitons are highly nonclassical, quantum correlated states 
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consisting of an even number of gap partides. This suggests the remarkable possibility that 
a bosonic frequency gap medium, while impervious to classical linear wave propagation may 
allow propagation of certain correlated quantum excitations. Multi-particle gap solitons may 
be generated by both nonlinearly exciting an impurity atom and Dicke superradiance from 
a collection of these excited àtoms. Unlike the single excitation which can tunnel a distance 
given by the classical penetration length within the gap, the paired excitations as well as the 
resulting heavy gap solitons can propagate freely through the gap of this harmònic medium. 
The Bethe ansatz solution which we have presented relied on the existence of an isotropic 
polariton dispersion relation. In a real PBG material, the photon dispersion relation is 
highly anisotropic. The propagation of quantum information within a PBG material, in this 
manner, would be of considerable importance in such applications as quantum computing 
\ïï\. 
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